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We develop a quantitative analytic theory that accurately describes the odd-even effect observed 
experimentally in a one-dimensional, trapped Fermi gas with a small number of particles [G. Ziirn 
et al., Phys. Rev. Lett. Ill, 175302 (2013)]. We find that the underlying physics is similar to 
the parity effect known to exist in ultrasmall mesoscopic superconducting grains and atomic nuclei. 

However, in contrast to superconducting nanograins, the density (Hartree) correction dominates over 
the superconducting pairing fluctuations and leads to a much more pronounced odd-even effect in the 
mesoscopic, trapped Fermi gas. We calculate the corresponding parity parameter and separation 
energy using both perturbation theory and a path integral framework in the mesoscopic limit, 
generalized to account for the effects of the trap, pairing fluctuations, and Hartree corrections. Our 
results are in an excellent quantitative agreement with experimental data and exact diagonalization. 

Finally, we discuss a few-to-many particle crossover between the perturbative mesoscopic regime and 
non-perturbative many-body physics that the system approaches in the thermodynamic limit. 

PACS numbers: 67.85.Lm,68.65.-k 


Our understanding of quantum systems is usually 
firmly rooted in either a few-body picture, where exact 
solutions of a few-particle Schrodinger equation exist, or 
a many-body picture, where the system can be described 
in a (quantum) statistical framework. In between these 
two limits lies the mesoscopic regime, where finite par¬ 
ticle number and confinement have a strong effect on 
the system’s properties. Mesoscopic systems occur nat¬ 
urally, for example, in nuclear physics, where a finite 
number of protons and neutrons form an atomic nucleus, 
or they can be engineered, such as in semiconducting 
quantum dots mm or superconducting nanograins El- 
E]. For attractively interacting mesoscopic Fermi sys¬ 
tems, a key effect is that the ground-state energy is not 
a strictly convex function of the particle number, but 
the interaction can cause some configurations to have 
lower binding energy (and thus enhanced stability) rel¬ 
ative to others mm- For example, this implies an en¬ 
hanced stability of nuclei with a “magic number” of con¬ 
stituents. A related effect exists for superconducting 
nanograins [ziiaiini: the binding energy of systems with 
an even number of spin-up and spin-down fermions (even- 
number parity) is enhanced compared to odd particle 
number systems with an unpaired fermion (odd-number 
parity). This parity effect is a hallmark of mesoscopic 
superconductor systems and can be quantified by the so- 
called parity or “Matveev-Larkin” parameter [ziisiiin], 
which denotes the excess energy of an odd parity state 
relative to the mean of the neighboring fully paired even 
parity states: 

Ap = E21+1 — - {E21 + E21+2 )) ( 1 ) 

where E 21+1 denotes the ground-state energy of a fermion 
system with odd total particle number 2^ -|- 1. For non¬ 
interacting systems, the parity parameter 0 vanishes, 


and it is positive if there is a parity effect. 

In this Rapid Communication, we study mesoscopic 
one-dimensional Fermi quantum gases, and establish a 
rigorous connection with well-known mesoscopic super¬ 
conducting systems. Our work is motivated by recent 
progress in quantum gas experiments which can deter¬ 
ministically prepare systems of few fermions in a har¬ 
monic one-dimensional trap m- These systems were 
studied for repulsive [T^] and attractive m interactions 
and spin-balanced m and spin-imbalanced configura¬ 
tions [14j , and used to simulate models of quantum mag¬ 
netism dSHn]. Motivated by a recent experiment [13) . 
here we study a spin-balanced few-fermion system with 
attractive interaction in a harmonic trap, i.e., ensembles 
which contain an equal number of spin-up and spin-down 
fermions for a total even particle number, and a single 
unpaired fermion for an odd total particle number. In 
Ref. [ini, following the preparation of an ensemble with 
a definite particle number, the trapping potential was 
tilted for a variable time, allowing fermions to tunnel out 
of the trap. From the tunneling times obtained in the ex¬ 
periment, a separation energy was extracted HaiSHlD], 
which is defined as 

E,,piN) = [En - E%) - [En-i - <_i), (2) 

where E^ {E%) is the ground-state energy of the inter¬ 
acting (noninteracting) system with N particles. At zero 
temperature, E^ is obtained by filling the lowest har¬ 
monic oscillator levels up to the Fermi level: for even 
total particle number 21, the states j = 0,... ,l — 1 are 
occupied by pairs of spin-up and spin-down fermions. For 
odd total particle number 21 1, the level I contains an 

additional unpaired fermion. The parity effect is mani¬ 
fested in the separation energy in the form of an odd-even 
oscillation, where the separation energy of an odd par- 
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tide number state is smaller than the separation energy 
of an even particle number state. The experiment [13] 
has been analyzed theoretically using exact diagonaliza- 
tion for small particle number However, for 

larger numbers of particles, exact diagonalization is be¬ 
yond computational reach and different theoretical ap¬ 
proaches are necessary. Recent numerical works compute 
ground-state properties using Monte Carlo methods for 
even fermion numbers up to = 20 |24| and coupled- 
cluster methods for up to = 80 [25l [26] . In this pa¬ 
per, we employ analytical methods, which allow a direct 
physical interpretations of the experimental results and 
provide complementary information to numerical works. 
Pairing in higher dimensions has been considered in |27b 

l 3 Dj. 


of a large system size L and large particle number N with 
constant density n = N/L, the parity parameter corre¬ 
sponds to half the spin gap, which for small interaction 
strength is Ap = f [31], where 
7hom = 'mgi /u <C 1 is the interaction strength of the ho¬ 
mogeneous system. For the trapped system, we expect 
that in the macroscopic limit of large particle number, 
the parity parameter is (in the Thomas-Fermi approxi¬ 
mation) given by its minimum value at the trap center 
where the local density is Uq = 2y/N /Trlho- This gives a 
parity parameter [25] 

A = Ap(A^oo) = —(5) 

TT V 


In the following, we analyze the mesoscopic pairing 
problem, focusing on the weak-interaction limit which 
corresponds to the experimental situation [13] . The par¬ 
ity parameter takes a fundamentally distinct form in the 
few-body and the many-body limits, interpolating be¬ 
tween a simple perturbative form and a manifestly non- 
perturbative many-body expression. We estimate a crit¬ 
ical particle number which marks the crossover between 
the mesoscopic and the macroscopic regime, finding that 
this quantity scales exponentially with the interaction 
strength, which suggests that the mesoscopic description 
persists over a wide range of particle number. Our the¬ 
ory is in accurate quantitative agreement with the ex¬ 
periment m and provides a theoretical framework to 
study the mesoscopic regime where A 3> 1 , which is of 
fundamental interest to understand the emergence of su¬ 
perfluidity and superconductivity in physical systems. 

The Hamiltonian of a two-component Fermi gas in one 
dimension is (we set H = 1) 


H = 






(3) 

Here, ipaix) annihilates a fermion at x with mass m and 
spin cr, V{x) = vtuxp'x^ 12 is the harmonic trapping poten¬ 
tial with frequency w, and > 0 is related to the effec¬ 
tive attractive scattering length oi via gi = 2lma\. We 
write the continuum model ([^ in an oscillator basis by 
expanding the fermion operators in terms of simple har¬ 
monic oscillator states '"^tiere 

4>j{x) is a normalized harmonic oscillator wavefunction 
with energy Sj = uj{j + 1 / 2 ) and the operator Cja- annihi¬ 
lates a fermion with spin a in state j. The Hamiltonian 
in oscillator space is 


ja ijkl 


where Zho = ^/ijrnuj denotes the harmonic oscillator 
length. The coupling is now state-dependent with an 
effective interaction strength set by the overlap integral 
’^ijkl — Zho J* dx • 

The theory in Eq. (1^ can be solved in the absence of a 
trapping potential [341133] . In the thermodynamic limit 


where the dimensionless interaction strength is 

'^91 

7 = — 7= -• (6) 

2ViVa;Zho 

Equation ([^ is a manifestly nonperturbative expression. 
Note that despite the exponential suppression with 7 , the 
parity parameter Ap scales with the Fermi energy. The 
macroscopic limit is therefore characterized by A ^ w. 
By contrast, in the mesoscopic limit of small particle 
number where A to, we expect simple perturbation 
theory to hold. This is reminiscent of the Anderson cri¬ 
terion that marks the vanishing of superconductivity if 
the level spacing of a grain is larger than the bulk super¬ 
conducting gap [34] . Clearly, the crossover from a few to 
many particles is manifest in the parity parameter. The 
expression ([^ is extensive with particle number for con¬ 
stant 7 , indicating that the crossover should be studied 
while keeping 7 fixed, i.e., imposing gi ~ '/N. In the 
following, we consider the regime where 7 ^ 1 . 

We proceed by analyzing the theory of Eq. 0 in the 
weak-interaction limit to first order in the coupling g, ap¬ 
plicable to the mesoscopic regime A <C w. To this leading 
order, the ground-state energy is given by the expecta¬ 
tion value of Eq. 0 with respect to the noninteracting 
ground state. The separation energy in Eq. (|^ is thus 


i-i 


A'sep(2Z -f 1) — 5lZjjo 'y ( Wji 

(7) 

3=0 


1-1 


Esepi‘^1) = ~9l^ho 

(8) 

3=0 



which corresponds to the interaction energy of a sin¬ 
gle fermion in the outermost level interacting with the 
fermions in the lower shells. In Eqs. Q and ([^, we 
define the diagonal coupling Wij = Wijji, which can be 
determined in closed analytical form [SB]. Note that the 
perturbative interaction correction is due to a mean-field 
shift of the single-particle energies. In Eig.[^ we show the 
results of Eqs. 0 and (|^ for the separation energy along 
with the experimental measurement [13] (black error 
bars) for an interaction strength gi/wZho = 0.45, which 
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FIG. 1. (color online) (a) Separation energy as a function of 
particle number N for interaction strength gi/lhoOJ = 0.45. 
The leading-order perturbative solution of Eqs. 0 and Q 
are indicated by blue circles and the path integral fluctuation 
correction [Eq. ( |25[ )] by red diamonds. Lines are a guide for 
the eye. Our result is in excellent agreement with the exper¬ 
imental results by Ziirn et al. [13], which are shown as green 
error bars. In addition, we compare with results from exact 
diagonalization m, indicated by gray bars, (b) Orange line: 
separation energy including only pairing interactions. The 
Hartree density interaction is essential to fit the experimental 
data. 


corresponds to the value used in the experiment |I3j . Re¬ 
markably, the perturbative result provides a very accu¬ 
rate description of the experimental data and is also in 
very good agreement with results from a numerical exact 
diagonalization of the Hamiltonian (|4| (gray bars) jUj. 
The parity parameter is 


Ap(A < w) 


giwii r [u 

—uj ~ A, / —, 
2cj/}^o V A 


(9) 


where A was defined in Eq. (§. 

To gain insight into the physical mechanisms con¬ 
tributing to the separation energy and the parity param¬ 
eter, we assume that for weak interactions, pairing takes 
place predominantly within a harmonic oscillator shell, 
i.e., that the ground-state properties can be described as 
excitations of paired levels: levels are either occupied by 
a pair of fermions or empty. This implies that only the 
interaction terms that connect two levels are retained: 
Wij = Wiijj = Wijij = Wijji |3S|. The effective Hamilto¬ 


nian takes the form: 

Hes= Y 

n<7 ij 

-5iC Y 

There are two interaction terms: The first, which we call 
the pairing term, destroys a pair of spin-up and spin- 
down fermions in one oscillator level and creates a pair 
in a different one. The second, which we refer to as the 
Hartree term, does not create excitations but provides 
a density-dependent energy shift to the single-particle 
levels (note that the perturbative result is due to this 
type of interaction). A third possible interaction term 
which exchanges the spin between two simply occupied 
levels {Fock term) does not contribute to the balanced 
system that we consider. Compared to pairing models 
used for superconducting nanograins, the pairing interac¬ 
tion takes a more complicated level-dependent form and 
involves an additional Hartree term, which is in fact es¬ 
sential to describing the experimental data of Ref. m- 
Figure [^b) shows the leading order prediction for the 
separation energy only taking into account the pairing in¬ 
teraction. As is apparent from the figure, this prediction 
is in complete disagreement with the experiment. Note 
that the while the Hartree term is crucial for a correct 
description of the separation energy, it does not affect 
the parity parameter of Eq. ([^ to leading order. 

We obtain the ground-state energy for fixed particle 
number from the limit [7] 

lim Heff = min(E 7 v — fJ-N), (11) 

/3—>-oo N 


where is the free energy of the system obtained from 
the grand canonical partition function Z = = 

Xj. q’j. dgnotes the trace over all 

many-particle eigenstates), i.e., the grand canonical en¬ 
semble projects onto a sector with definite particle num¬ 
ber. However, because of the parity effect of Eq. ([^, the 
prescription © only allows us to study configurations 
with even particle number. Nevertheless, we can relate 
the ground-state energy of a system with an odd number 
2/ -I- 1 of fermions to the ground-state energy of a sys¬ 
tem with an even number 21: since the Hamiltonian (10) 


only couples fully occupied or empty levels, the unpaired 
orbital of an odd-particle number state does not partic¬ 
ipate in the interaction and decouples; i.e., it effectively 
blocks a level from the Hilbert space. Hence jiiiaiinj 


E2I + I — Sl+l + E21, 


( 12 ) 


where the energy E 21 is computed without the blocked 
level 1. To analyze the effective theory of Eq. (10), we 


eliminate the quartic interaction terms using a double 
Hubbard-Stratonovich transformation in both the den¬ 
sity and pairing channel, which introduces three auxil¬ 
iary fields Ai, A*, and AT^. To this end, we define the 
operators q° = qt = + 
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and ( 7 j = f(cLc-| — CjiCj^). The Hamiltonian (10) takes 


‘"it 4 WtWt 


T = Y.(>i 

3 

Va = gi^ho ^ 


i_ — 14- — Vo, where 


+ Y 

i 

(13) 

Wijqfqf, 

(14) 


with = Ej—^. The last constant term in Eq. ( |13[ ) arises 


from a fermion commutator. Since the symmetric ma¬ 
trix {w~^)ij/gil^^ is positive definite, the four-fermion 
interaction terms can be reduced using three standard 
Hubbard-Stratonovich transformations m for the op¬ 
erators gf introducing conjugate real fields xf, where 
a = 0,±. Identifying Kj = Xj and Aj = x^ + ixj, the 
partition function reads: 


Z = 


X exp 


1 

37 


l[VA3T)VA*iT)VK3r) 


Tr[t/A(/3,0)] 


-pc- dr 

Jo 


(,W 


^3 


9il 


^ (A*A, + K,K,) 


ho 


(15) 


where C = + Ci) with - Kj, Af is the 

path integral normalization 




X exp 


Y[VA3T)VAUr)VK3r) 

■ r,i 

- f dr Y, (A* A, + K,Kj) 


, (16) 


and 

Ua{P, 0 ) = Tt exp 


- drYx, 

^0 


^3 


(17) 

with X3 = iCjj,c]Y M- 

We first consider the saddle-point approximation and 


minimize the Euclidean action in Eq. (15) with respect 


to Ki and A^. To this end, we first integrate out the 
fermions in the partition function, which gives the effec¬ 
tive action 


guy 

ih.. ^ 


Aj 

^4 2Ei ’ 


j ~,j 2 t] ; where Ej = Tj + A^. The solution of 
these equations determines the value of the Hartree field 
Ki and the gap Aj at the saddle point. Note that for an 
even particle number, the saddle-point equations corre¬ 
spond to the solution of a BCS pairing ansatz [MiiMiiin] 
with coherence factors vf = l — uJ = [l — {^i — Ki)/Ei]/2. 

For small 7 , the only saddle-point solution for Aj 
corresponds to Aj = 0 , which implies vanishing off- 
diagonal long-range order and absence of superfluidity 
in the weak-coupling regime. This is the fluctuation- 
dominated regime where A defined in Eq. © is much 
smaller than the harmonic oscillator level spacing ui [3]. 
The Hartree fields are given by Ki = 
which correspond to the single particle energy shift com¬ 
puted to leading order in perturbation theory in g using 
the noninteracting ground state of the Fermi gas. Using 
the identity in Eq. (11), the saddle-point contribution to 
the ground-state energy of an system with even N = 21 
particle number is E 21 = 2 X)j‘=o ^3 “ 9i Z)iyLo + 
^ 'Tli'^ii- Interestingly, this is not equal to the pertur¬ 
bative ground-state energy. The last term arises from 
the commutator term in Eq. (13). Despite the saddle 
point Ai being zero, fluctuations of the pairing field can 
make an important contribution, and they are computed 
in the remainder of this paper. It turns out that these 
pairing fluctuations contain a 0{gi) correction that can¬ 
cels the last term in the saddle-point contribution to E 21 
and reproduces the perturbative result. 

To consider the effect of fluctuations around the 
saddle-point solution. We write Ki ^ Ki + SKi and 
Ai SAi, and expand the action in Eq. (15) to sec¬ 


ond order in 6Ki and SAi. It turns out that there 
is no correction due to fluctuations of the Hartree 
fields SKi. The partition function can be written as 
Z = ZgpZsA, with Zgp the saddle-point contribution 
and ZsA = [detr]“^ the resulting quadratic functional 
integral in SAi, which can be exactly evaluated in 
terms of the functional determinant [3 [37]. To de¬ 
rive this result, we expand the perturbation in Mat- 
subara space SK^t) = YJ2iuj,^^~^‘^"'^SKi{iuJn), and 
(5Ai(r) = -^J2iujr, e“*‘^"’’dAi(za;„). The effective action 


is 


5eff[{A„X,}] = -)_^trln[-Go-j] 


+ dr 

Jo 


^3 




93 


-K,K, 


ho 




9i 


2 ^ 


Wii 


where the trace runs over Matsubara indices, 
matrix element of Gqj is given by (*a;„|Gojjia;„/) = 

Sn,n'CQj3uJn) — ^n,n' “ Sn,n'' 



-S'eff [{6'i}] 


fp 


+ dr 

The 

Jo 


Varying the action S'eff in Eq. (18) with respect to Kj, 
Aj, and /r, we obtain the mean-field saddle-point solu¬ 
tion defined by Ki = YJj'^ij (l ~ = 


(w l)ij 
-1 

ho 


93 


SA*SAj + {K, + SKJ) 


x{K,+SK,)\ 


(19) 


where we separate a fluctuation part 14 from the Green?s 
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function Goj - The matrix element of Vj is given by 


Using 


{iLOn\Vj\iuJn') = Vj{iuJn - iu)n’) 

= - \= 6Kj{i0Jn - iUJn')cr3 + SAj{i0Jn - iu}n')cr~'~ 

VP L 


+SA*{iu}n — iujn>)a- 


( 20 ) 


LXJ ^ 

- = EE ( 21 ) 

3 3 "=0 


we expand the effective action (191 to second order in Vj. 


The functional integral in terms of 6Kj and SAj can then 
be performed exactly. The partition function involving 


SKj reads: 


ZsK = 


Nk J 

X exp 


Y[vSKk 

k -■ 

- EE = 0 )f 2 ^ 


,_i -‘-V --- It / - j \ - / V ^ - ^ ,_i --J ■ 1^1 


0 


- 1 


.( 22 ) 


However, the zero-frequency contribution (second term) vanishes since Kj is given by Kj = ■ The 

remaining quadratic term is irrelevant since it does not involve any single-particle energies. The partition function 
involving 5Aj reads (discarding an irrelevant constant term) [71 137]: 




\{V5AIV5A, 


= det 


k 

-1 


exp 


lOJn 


X I 7-1 

Oij + QiIy^o 


sgn^, 

*? • 


-EE“;(- 

sinh f3^j 


lUJn IJ 


n 


sinh ‘ 


9il 


ho 


sgnQ 

1 * O G 


5Aj{iu3n) 


(23) 


r 


where by { 2 ^^} we denote the eigenvalues of Aij = 
‘^S. 3 ^i 3 - 9 iWij sgn^j, i.e., det( 2 |fc/ - H) = 0 , or, respec¬ 
tively, 


detf (5: 


9ikl sgn^j 


= 1-b 


2 

9iK^ 


‘ho 


6 - ^3 

sgn^i 

ik - Ci 


E- 


= 0 , 


(24) 


where the second term holds for small corrections. Writ¬ 
ing ^j = + 6 ^j and expanding in S^j , the fluctuation 

correction to the free energy at zero temperature is: 


6n = -EnZiA = E'^Osgn^j, 


(25) 


where 
(5^2 = 


giwu sgn^, 

2 ^ho 


1 - 


9l WjjSgR^j 

2^ho E 0 


-1 


(26) 


Hence, SE 21 = J2j ^^3 sgn^i and SE 21+1 
^^7 sgn^', where and 6 ^' are computed 
in Eq. (26) with the ^-th level excluded. 


as 


The fluctuation correction to the separation energy 
and the parity parameter can be read off directly from the 
definitions 0 and (H- Note that the combined saddle- 
point and fluctuation correction contains the leading or¬ 
der perturbative result (see Fig. where the separation 


energy is indicated by the red dashed line and the dia¬ 
mond symbol). There is a small quantitative correction 
which improves the agreement with the exact diagonal- 
ization results by D’Amico et al. |7T]. The fluctuation 
correction (26) is similar to the one encountered in su¬ 


perconducting nanograins in the limit where the super¬ 
conducting gap A is much smaller than the level spacing 

Se |7]. 

Interestingly, our analytical procedure also allows us 
to identify the critical particle number that marks the 
crossover between the few-body regime A VC w and the 
many-body regime A 3> w. The boundary of the meso¬ 
scopic regime is marked by a breakdown of the expan¬ 
sion (26). For large particle number, we can replace the 


harmonic oscillator matrix element by its semiclassical 
expression wjj ^ and convert the summation to an 

integral. This gives 


9\ 

2 ^ho 


E 


- ^2 


gi 

tolho 


Aln 22 V, 


7 In 2 A, (27) 


indicating that by successively increasing particle num¬ 
ber, the few-body expansion loses validity at a critical 
particle number Nc ~ In this case, the bulk par¬ 

ity parameter ([^ is comparable to the level spacing w, 
which is a corresponding criterion as in superconducting 
nanograins |^. Note that while the few-body to many- 
body crossover is manifested in the parity parameter at 
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leading order, the ground-state energy is dominated by 
a Hartree mean-field term and will be less sensitive to 
the crossover. From the perspective of superconducting 
nanograin systems, such a predominance of the Hartree 
contribution over the fluctuation correction is an unex¬ 
pected effect [HISlIlo]. Therefore, our findings prompt 
a revision of both the theoretical modeling of nanograins 
and the related experimental results EH- 

While the BCS pairing model can be solved ex¬ 
actly [4lVl48] , this is not the case for the model Q or the 
reduced pairing Hamiltonian (101. However, it would be 
interesting to explore if the theory could be approximated 
by a generalized Richardson-Gaudin model [iH] , 

In summary, we have computed the ground-state en¬ 
ergy, the separation energy, and the parity parameter for 
a trapped one-dimensional Fermi gas with weak attrac¬ 
tive interaction. We have used an insightful path-integral 
formalism which allows us to make useful connections 
with other physical systems (i.e., mesoscopic supercon¬ 
ductors). The parity parameter serves as an order pa¬ 


rameter that displays a fundamentally distinct behavior 
in the mesoscopic and macroscopic regimes, and we es¬ 
tablish that the mesoscopic description persists for a wide 
range of particle number. Our results provide a quanti¬ 
tative description of the recent experiment m- A path 
integral treatment indicates that the ground-state energy 
and the parity effect are dominated by a Hartree mean 
field contribntion, with BCS pairing flnctuations provid¬ 
ing a subleading correction to this result. 
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